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1. Introduction 

The integral quantum Hall resistance [1] 
R H (i)=V H (i)/I x =h/(eH) is observed when the 
longitudinal resistance R x = V x /I x of the two-di- 
mensional electron gas is very small. Here F H (0 is 
the Hall voltage of the / th plateau, / is an integer, 
and I x is the current through the sample. There is a 
critical current 7 C above which R x rapidly increases 
by several orders of magnitude [2,3]. R x becomes 
finite as one approaches the critical current. This is 
referred to as breakdown. 

We present results of a breakdown experiment in 
which sets of discrete V x signals are observed. We 
propose an explanation based on an extension of 
the quasi-elastic inter-Landau level scattering 
(QUILLS) mechanism of Eaves and Sheard [4]. 



2. The Experiment 

Our samples [3] were GaAs/Al^Ga^As het- 
erostructures grown by molecular beam epitaxy 
with x=0.29. The samples are designated as 



GaAs(7) and GaAs{8). They have zero magnetic 
field mobilities of -100,000 cmV(V-s) at 4.2 K. 
They exhibit excellent integral quantum Hall effect 
properties. The inset of figure 1 shows the geome- 
try of the samples. They are 4.6 mm long and 0.4 
mm wide. Contact was made to the two-dimen- 
sional electron gas at points 1, 2, 3, and 4, as well as 
the source S and the drain D. 

An important feature of the GaAs(7) sample is 
that the breakdown occurs somewhere within the 
longitudinal shaded region between the middle 
Hall probe pair 3,4 and the outer Hall probe pair 
1,2 of figure 1, but not within the transverse region 
of either Hall probe pair. This crucial fact can be 
explained by examining figure 1 . The minimum V x 
signal, V? m , measured between probe pairs 2 and 4, 
increases by a factor of 10 7 between 7* =25 /xA and 
7* =3 70 jiiA, but the quantized Hall resistance R H 
for the /=4 plateau decreases by only one part in 
10 7 for probe pair 3,4 at 370 /xA, and by only six 
parts in 10 7 for probe pair 1,2. These changes in R H 
are only about 0.01% of what was expected from 
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Figure 1. Current dependence of ARh/Rh and VJ?'" for the i=A 
plateau of the GaAs(7) sample at 1.2 K. The inset shows the 
sample geometry. 

the mixing of V x into V H due to the known mis- 
alignment of the Hall probes. We therefore know 
the general region where the breakdown occurs, 
and that the / =4 Hall resistance is accurately quan- 
tized on both sides of the spacial breakdown re- 
gion. 

The critical current at which V? m starts to rise 
abruptly is 7 C = 340 juA. Figure 2 shows a V x vs B 
curve for GaAs(7) at 7^ = 300 ju,A, well below that 
critical current. This curve exhibits reproducible 
structures. Figure 3, (a) and (b), respectively, show 
digital oscilloscope displays of the time-depen- 
dence and time-averaged distributions of the values 
of V x obtained at points A, B, and C of figure 2. 
There are clearly distinct dc voltage levels and 
switching between levels. Figure 3 clearly shows 
that the character of these levels changes with 
magnetic field. Also, V x is in only one state at any 
given time. It remains in that state until electrical 
noise or other noise processes induces it to switch 
to another state. Similar switching between dc 
voltage levels occurs in the GaAs(8) sample. The 
V x curve of figure 2 masks this switching by dis- 
playing time-averaged values of these dc voltage 
levels. 

We make an initial attempt to interpret our ob- 
servations by using a modified version of the 
QUILLS model of Eaves and Sheard [4] because 
that is a model which seems to provide a satisfac- 
tory explanation. Eaves and Sheard used their 
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Figure 2. Time-averaged V x signals for the ;=4 plateau of 
GaAs(7) at 1.2 K and 300 jxA. The number of inter-Landau 
level transitions (N'—N) is indicated in brackets. 

model to interpret the experimental data of Bliek et 
al. [5] who made breakdown measurements on 
GaAs/AlGaAs samples in which Hall potential 
probe sets were placed on either side of a 1 ju,m 
wide constriction. Their V x vs B curves show fea- 
tures similar to our data, and their Hall voltages 
were also quantized. Our samples have no geomet- 
rical constriction; they are instead 400 ju,m wide. 
But we propose that the nearly dissipationless con- 
duction channel is very narrow in our samples as 
the critical current is approached, and that the con- 
duction channel becomes entirely pinched-off and 
dissipative above the critical current. Our samples 
therefore may have effective geometries similar to 
those of Bliek et al. [5] near breakdown. We assume 
the global current-carrying equipotentials of the 
percolation model [6], as strongly suggested by ex- 
periments which measured the potential distribu- 
tion across quantum Hall devices [7]. 



3. The QUILLS Model 

To explain our data we first extend the QUILLS 
model of Eaves and Sheard [4] to include transi- 
tions between non-contiguous Landau levels. Parti- 
cles of positive charge q and reduced mass m* 
(0.068 times the free electron mass in GaAs) move 
with velocity v x in the positive x direction through 
a constriction of average length L x and average 
width L y as shown in figure 4. The quantities L x 
and L y can be less than the length L and width Wof 
a geometrical constriction. The average electric 
field within the constriction and the average mag- 
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Figure 3. (a) Digital oscilloscope recordings of the time-dependent V x signals at the positions labeled A, 
B, and C in figure 2. (b) Digital oscilloscope recordings of the time-averaged distribution of the values 
of V x at the positions A, B, and C. 




Figure 4. Motion of a charged particle through a constriction. 
The dynamical length and width L x and L y may be smaller than 
the physical length and width L and W. 

netic field is E y and B zt respectively. Therefore, 
v x =E y /B 2 . The Hamiltonian for this system, ne- 
glecting spin-splitting and scattering, is H=(\/ 
2m*)(p—qA) 2 +qV y in SI units, where p is the 
momentum operator —ilN, A is the magnetic vec- 
tor potential and V y — —yE y . Using the Landau 
gauge A x = —yB z and A y =A 2 =0, one obtains the 
normalized eigenfuctions i|f and energy eigenvalues 
% to Schodinger's equation jfy=g?i|i: 



tyN,k x ( x >Z)—/T y/2 e ' ' x (2 N M\) U2 C'j7-V /4 e 2 Hn(H) 



V«*x(*»0=^i/2 e' M <M£> 



^;v(yo)= (N+ -)Ha c -qyoE„ +-m *v 2 x , 

where TV is the Landau level number, oi c =qB z /m* 
is the cyclotron angular frequency, f B =(fi/qB z ) i/2 
is the magnetic length, ya—(v x /o) c — <f B k x ) is the y 
coordinate of the center of motion of each cy- 
clotron orbital, k x = 2TrN k /L x is the x component 
of the wavevector of a cyclotron orbital state 
(N,Nk), N k is the integer quantum number of the 
wavevector of that state, %=(y— y a )/f B , and 
H N (t£) is a Hermite polynomial. 

Figure 5 shows the energy eigenvalues for 
Landau level TV as a function of y. The slope of the 
lines is —qE y . The eigenvalues are equally spaced, 
with separation Ay =2Trf 2 B /L x . The maximum 
number of allowed extended states (N,N k ) is (qBJ 
h^LJ-.y for level N. The total number of allowed 
extended states per unit area is n s =i(qB 2 /h). 
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Figure 5. Total energy eigenvalues 8 s N as a function of position 
y across the constriction for each Landau level N. Every eigen- 
value of level N has a unique quantum number A'*. The figure 
shows a QUILLS transition from level N to N' and an associ- 
ated acoustic phonon. The decay, either directly down to level 
N ox through an intermediate level N", and its associated opti- 
cal phonon, are also shown. 

If L y becomes small enough and E y large enough 
then it may be possible for a particle to make a 
transition from state (N,N k ) to an empty state 
(N',N k ) at a lower total energy, as shown in figure 
5. The particle moves across the sample in this 
model from position y to y{>. Energy and momen- 
tum must be conserved in the transition. Therefore, 
an acoustic phonon of wavevector 



K x =M(N k -M)= 



G>c 



(V x -V s ) 



(N'-N) 



(1) 



and energy ^ phonon =-fia s =fw s K x is emitted in the x 
direction, where v s is the velocity of sound in that 
medium (-2.47X 10 3 m/s in GaAs[8]). 

Only transitions in which the ratio (y£ — y )L x / 
(27r/|) is an integer number (N k —N k ) are al- 
lowed. One can obtain an equation for (N k —N' k ) 
by estimating the value of (yi —yo)- This is accom- 
plished by noting that the spatial extent of 
tyN,k x (x>%) or ^n(£) is approximated well by the am- 
plitude of motion of a classical harmonic oscillator: 
A N = ^ b^IN+X. Transitions can commence 
when the wavefunctions just begin to significantly 
overlap: 



(yi-y^(A N +A}t)=/ B ^ 



N.N 1 



(2) 



where gf r =V2N+ 1 +V22V + 1. Figure 6 shows 
such an overlap between <J> (§) and <J> 12 (£). Thus 
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Figure 6. The wavefunctions <I>o(€) and < t>i2(€)> shown overlap- 
ping at the classical harmonic oscillator separations § and £i 2 . 
The dashed lines represent the classical harmonic oscillator 
probability distributions. 
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There is another condition for the transition: 
(N'—N) must also be an integer. From conserva- 
tion of energy, the electric field E y = V H /L y is 



Ey = [(N' - JV)fc»C + fc» J/?0\) -J>0)- 



(4) 



The simultaneous integer conditions (N k —N k ) and 
(N'—N), given by eqs (3) and (4), occur only at 
particular values of the current and magnetic field. 
There may be many intervening Landau levels be- 
tween N and N' in which these conditions are not 
satisfied. Notice that the filled states of Landau 
level N are at a higher total energy than those of 
the unoccupied N' level. The large Hall electric 
field E y has induced a population inversion for 
QUILLS transitions. 

After the particles make the N to N' QUILLS 
transitions they then decay, by optical phonon 
emission [9], back down to the original Landau 
level N, either directly or by cascading through 
intermediate levels N", as shown in figure 5. The 
optical phonon transitions probably occur just out- 
side of the constriction. One can calculate the 
voltage signal V x from the electrical power dissi- 
pated, P=I X R X =I X V X , when the particles return to 
the lowest Landau level during the average time 
taken to traverse the constriction: 






(5) 
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where I x =qn & v x L y , and n'/n s is the fraction of con- 
ducting particles that make the transition. We as- 
sume both spin states can undergo transitions. 
Therefore iVis and «7« S <1 for the i=2 plateau. 
N is either or 1 and «7rc s <l/2 for the /=4 
plateau. If all the electrons make the transition, 
then n'/n s is 1 and 1/2 for the i=2 and 4 plateaus. 



4. Analysis of Our Data 

We next apply this QUILLS model to our V x vs 
B data at 300 jaA. Equation (5) is used first to ob- 
tain the values of (N'—N) by ignoring the small 
term fta s . The electric field E y is then calculated 
from eqs (4) and (2) for specific values of V x and B 2 . 
E y ranges between 1.6-3.9X10 6 V/m across the 
quantized Hall resistance plateau. The constriction 
width L y = Vn/Ey =I X R H /E y varies between 0.5-1.2 
fim across the plateau, and is narrower on either 
side of the V x zzO region. The range of v x =E y /B z is 
2.8-7.0X10 5 m/s across the plateau, and v x is 
~200\) s . Also, fua t =fivJK x , calculated from eq (1), 
is -0.5% of (N'-N)Ha c . The dissipation of the 
QUILLS transitions is therefore very small indeed. 
The current density J x =I x /Ly=E y /(h/Aq 2 ) varies 
between 240-600 A/m. This is an astonishingly 
large number compared with that found in [5], but 
the current has a better opportunity of finding a 
dissipationless path in wide high-quality samples. 
Also, qV u ^200Hoi c . This is very large compared 
with our deduced (N'—N)fio) c values, and is in an 
entirely different regime than that of Kirtley et al. 
[10] who found small integer values of qV H /fia> c . 
(N k —Nk) is the only quantity that can not be de- 
duced from the data because the length of the con- 
striction L x is unknown. If L x is assumed to be 
comparable to the width, i.e., about 1 jxm, then 
(jV^-iVOi^O from eq (3). Finally, eq (5) predicts 
that the quantity (n/n')(m*/-fi)V x /[(N'-N)B z \ 
should be nearly equal to 1. This quantity is within 
3% of 1 for the (N'—N) transitions shown in fig- 
ure 2, and is within our experimental accuracy. All 
of the electrons in Landau level N apparently make 
the transition to level N'. 

Multiple values of V x sometimes occur at certain 
values of B z , as for example at points B and C of 
figures 2 and 3. In such cases the above equations 
yield a different value of L y for each value of V x . A 
more realistic approach in those situations is to use 
the smallest value of L y obtained from the largest 
value of V x for that magnetic field. The values of L y 
and E y are then constant for all QUILLS transi- 
tions observed at a given B 2 and I x . With these as- 



sumptions, the values of £, N ' M for the smaller V x 
tr ansition s ar e then less than the classical values 
V2AT+ 1 + V2AT + 1. This presents no problem be- 
cause the overlap integral between <&#(£) and 
<1>ap(£') becomes somewhat larger inside of the clas- 
sical value and then falls off slowly over a wide 
region as £ wv is reduced. 

The features labeled with the (N'—N) values 12, 
13, 15, 16, and 18 on the right hand side of the V x 
curve of figure 2 have very stable single-valued sig- 
nals. This is not always the case; for example, we 
see switching about the (N'—N)=\2 transition at 
position C and about the (N'—N)=2\ transition at 
position B. We believe that the switching is noise- 
induced. 

We clearly observe discrete levels of V x in figure 
3. It is difficult to understand why the overlap of 
the rather broad wavefunctions and the integer re- 
quirements of the apparently large values of 
(N'—N) and (N k —N' k ) would by themselves lead 
to preferential N to N' inter-Landau level transi- 
tions. An additional, unaccounted selection mecha- 
nism may be present. 

One has to be careful about the definition of crit- 
ical current. For example, in our experiment V x is 
large and there appears to always be dissipation at 
points B and C in figure 2. However, when looking 
at the discrete voltage states of those points in fig- 
ure 3 we see that for one of those states (V x ^0) the 
dissipation is negligible. Therefore the critical cur- 
rent has not yet been exceeded. 



5. Discussion 

Our interpretation of the QUILLS model for 
breakdown is consistent with the experiment of 
Bliek et al. [5], who made physical constrictions of 
order 1 jum width and lOjum length. They observed 
step-like structures in their V x vs B curves and em- 
pirically found equations for the quantized values 
of R x at those steps and for the magnetic field val- 
ues at which the structures occur. The empirical 
number n L —L x (2qB z /h) in of Bliek et al. [5] corre- 
sponds to n L =(N k —N! c )(<l7ry /2 /£cf in our formu- 
lation, whereas Eaves and Sheard [4] predicted 
that n L =(N k —N' k ). The equation R x =n R (n L / 
n )(h/2e 2 ) of [5] corresponds to our interpretation 
of QUILLS if n R =(n'/n s )gf/(ir) m and 
n =n s L x L y . Eaves and Sheard [4] predict that n R is 
an integer number. Neither n L nor n R are integers in 
our formulation. 

In our interpretation of the data of Bliek et al. [5] 
the values of (N'—N) vary between 3 and 13 for 
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their i = 1 plateau assuming, as they did, that their 
V x signal is not a time average of several Landau 
level transitions. They obtained values of n R be- 
tween 2 and 19. We find that the effective constric- 
tion of their sample L y varies between about 0.2-0.5 
jam, significantly less than the 1 jam physical width. 
If L x equals the physical 10.3 jam length L, then 
(N k —Nk) ranges between about 650-990. L x is 
probably much less than L , yielding smaller values 
of (N k —Nk). Bliek et al. [5] found that n L varied 
between 528 and 607. 

Sachrajda et al. [11] also used samples with a 
narrow constriction. They observed structures 
having the same triangular behavior as that of the 
critical current versus magnetic field plot of Bliek 
et al. [12]. They believed that this behavior did not 
agree with the QUILLS model. But we find from 
our data that the conduction channel narrows on 
either side of V x ~0, and because less current is 
required in a narrower channel to obtain the same 
electric field, our results show that the QUILLS 
model is consistent with that behavior. 

It has been proposed [13,14] that breakdown is 
due to emission of phonons in a manner analogous 
to the Cherenkov effect if v x >v s . In the QUILLS 
interpretation of our data o^^OOu, and the 
Cherenkov angle c = cos _1 (v s /'U c ) varies between 
89.5° and 89.8°. The dissipative voltage V x can be 
very small even when the velocity v x is apparently 
quite large. The experiment of Bliek et al. [5] pro- 
vides excellent proof that this may indeed be the 
case. Figure 1 of their paper shows that V x has 
structures on the sides of the i = 2 plateau and goes 
to zero in the center of the plateau for I x =45 juA, 
L y ^\ jam, and 2? r ;=6.4 T. Therefore, v x =E y / 
B z = (I x R H )/(L y B z )^26.7v s , and yet V x still goes to 
zero. 

We have considered the acoustic phonon emis- 
sion to be collimated in the x direction. However, 
the electron gas has a probability distribution in the 
z direction that extends over an interval Az which 
is about 50 A in GaAs/AlGaAs. Therefore the 
acoustic phonons may have K z ~ 1/Az components 
which are about 35% of the magnitude of the K x 
component in our experiment for the (N'—N)= 12 
transition. Phonon emission in the x-z plane would 
not have a significant effect on the results of this 
paper, but it might preclude the possibility of spon- 
taneously emitted phonons stimulating the other 
phonons to make population-inverted QUILLS 
transitions in a manner analogous to photon emis- 
sion of lasers. 



6. Conclusions 

We have observed spatially localized breakdown 
of the nearly dissipationless quantum Hall effect 
into a discrete set of dissipative states, and have 
interpreted these observations in terms of the 
QUILLS model. Theoretical analyses of QUILLS 
transition rates, similar to [13], and optical phonon 
decay rates would be very useful in determining 
the validity of QUILLS. Direct observation of the 
acoustic phonons would confirm the QUILLS pro- 
cess. There are other interesting features that could 
be studied: the decay process; the quantum heater 
property of R x ; and the fact that two crossed resis- 
tances, Rh and R x , are both quantized. We suggest 
that QUILLS transitions may be the dominant 
breakdown mechanism in high-quality integral 
quantum Hall effect samples, and that the discrete 
V x signals are indicators of those transitions. 
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